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The Kibble-Zurek mechanism demands an initial adiabatic stage before an impulse stage to have 
a frozen correlation length that generates topological defects in a cooling phase transition. Here we 
study such a driven critical dynamics but with an initial condition that is near the critical point 
and that is far away from equilibrium. In this case, there is no initial adiabatic stage at all and thus 
adiabaticity is broken. However, we show that there again exists a finite length scale arising from 
the driving that divides the evolution into three stages. A relaxation-finite-time scaling-adiabatic 
scenario is then proposed in place of the adiabatic-impulse-adiabatic scenario of the original Kibble- 
Zurek mechanism. A unified scaling theory, which combines finite-time scaling with critical initial 
slip, is developed to describe the universal behavior and is confirmed with numerical simulations of 
a two-dimensional classical Ising model. 
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The Kibble-Zurek mechanism (KZM)^“^ describes 
topological defect formation in driven critical dynam¬ 
ics in a variety of systems, ranging from classical^^® to 
quantum phase transitions^®'''^*’. Kibble first proposed it 
in cosmology’’® by identifying a frozen correlation length 
^ that renders spatially distant regions causally indepen¬ 
dent during the cooling of the universe from the big bang. 
Then, Zurek brought this proposal to condensed mat¬ 
ter physics and offered a method to compute the density 
of defects formed'^’’. As a system cannot always follow 
adiabatically the cooling of a finite rate R due to crit¬ 
ical slowing down near the critical point, its evolution 
from a temperature Tq, sufficiently higher than the criti¬ 
cal temperature Tc, can be divided into three sequential 
stages, an initial adiabatic stage, an impulse stage, and 
a final adiabatic stage below Tc- In the initial adiabatic 
regime, the correlation length ^ and the correlation time 
Qs grow as \£\~'' and respectively, as the distance to 
the critical point, e = T — Tc, is reduced, where v and 2 ; 
are the correlation-length and the dynamic critical expo¬ 
nents, respectively”. The boundaries between the stages 
are then determined by the frozen instant t at which the 
time interval before the transition, tc — t = e/R, equals 
Cs®’’, where tc = Sq/R is the time at e = 0. This leads 
to tc — i ~ 3,4^ where rp = z + l/i' is a rate 

exponent'^®. Upon assuming evolutionless in the middle 
impulse stage, i then determines ^ ~ and thus 

the defect density n ~ the KZ scaling’h^. 

Crucial in the derivation is the existence of the ini¬ 
tial adiabatic stage that gives rise to It results from 
the large eo = To — Tc and thus small (s- By contrast, 
whether the initial state is equilibrium or not is irrelevant 
as the system can quickly equilibrate once Cs is small. 
This has been confirmed by a lot of experiments and nu¬ 
merical simulations®'^'®®’'’®''^^. On the other hand, when 
£0 = 0 and the initial state is the equilibrium state there, 
it has been shown by an adiabatic perturbation method 
that the scaling of topological defects is consistent with 


the KZ scaling'^®’"’®’'’®. Yet, it is difficult to obtain the 
equilibrium state near the critical point due to critical 
slowing down. 

However, if £q is small and the initial state is not the 
equilibrium state at Eq, the equilibration of the system 
has to take a long time as the relaxation time is now 
macroscopically large. In this case, there is no initial 
adiabatic stage at all and thus adiabaticity is broken. 
Questions then arise as to whether there is still a ^ that 
generates topological defects, or, whether the KZM is 
still valid or not. Does universal behavior exist in this 
driving critical system with a nonequilibrium initial con¬ 
dition? If the answers are yes, then how is ^ determined 
and how does one describe the universal behavior, as the 
adiabatic-impulse-adiabatic scenario of the KZM cannot 
apparently be applied to this case? 

Relaxation of a nonequilibrium initial state near Eq is 
not a strange situation^®'®®. A well-known case is the 
critical initial slip®®’®'*, which was found in classical’®’®* 
and recently in quantum critical phenomena in imagi¬ 
nary time”®. When a system is quenched rapidly from 
a high temperature disordered state to near its critical 
point and relaxes, it has been found that the order pa¬ 
rameter M grows as M ^ M^t^ right after a microscopic 
timescale, where 9 is an independent initial-slip exponent 
and Mq is a small initial order parameter, which may be 
generated by an external field®®’®*’®®. As the initial state 
is derived from the disordered phase, it possesses only 
short-ranged correlations. However, finite-ranged corre¬ 
lations are irrelevant in the renormalization-group (RG) 
sense®®’®*. So, the initial state may be an equilibrium 
state of a Hamiltonian different from the system’s. 

A system that is driven by an external field including 
the temperature with a constant time rate R through its 
critical point is well described by the theory of finite¬ 
time scaling (FTS)®^’®®. It is a temporal counterpart of 
the well-know finite-size scaling*’ and is derived from the 
RG theory”’"’®®. FTS shows that there is a hnite timescale 
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Cd induced by an external driving; and when 

Cd is shorter than it dominates the evolution in an 
FTS regime. This indicates that the impulse regime of 
the KZM is just an FTS regime. Indeed, Cd is just i 
because at this instant Cd = Cs uud C is just the length 
scale corresponding to Q. Moreover, the scaling behavior 
in the evolutionless impulse regime and the KZ scaling 
are well described by In FTS, however, the 

initial state is, similar to the KZM, far away from the 
critical point and has thus no effects. 

Here, in order to describe the scaling behavior of a 
driven critical system with a nonequilibrium initial con¬ 
dition, we combine FTS with the critical initial slip. We 
shall show that there again exists in this case the finite 
timescale Cd and thus C- As a result, the KZM for topo¬ 
logical defect formation is still valid though adiabaticity 
is broken. However, its adiabatic-impulse-adiabatic sce¬ 
nario is now changed to a relaxation-FTS-adiabatic sce¬ 
nario, in which a nonequilibrium nonadiabatic relaxation 
stage replaces the original initial adiabatic stage. In this 
relaxation stage, the growing correlation time Q, which 
is different from Cs ^ dominates the evolution 

and Cd is subsidiary. Once Ci gets longer than Cd, the 
latter takes over and the system enters the FTS stage. 
This is the impulse stage of the KZM. However, in the 
KZ sense, both the relaxation and the FTS stages are 
impulse as both are nonadiabatic albeit due to different 
reasons, viz., the former arises from the initial conditions 
whereas the latter from the driving. When the system 
is driven to so far away from Tc that Cs becomes shorter 
than Cd, it crossovers into the adiabatic stage. 

As an appreciation of the results, we plot in Fig. 1 
the evolution of M for three different sets of the initial 
conditions. One sees that the evolution starting with a 
large |£o| and at Tc show qualitatively distinct behav¬ 
ior. Driving {R ^ 0) makes no appreciable difference 
at short times when the system starts with an uncorre¬ 
lated nonequilibrium initial state near to its Tc- In this 
stage, since the correlation length grows as Ci ~ and 
C, ~ ~ As Ci is shorter than Cd in short 

times, it dominates the dynamics and the stage is thus 
relaxational similar to the critical initial slip, while the 
external driving in only a perturbation. Note that this 
relaxation stage has nothing to do with the free relax¬ 
ation regime^^’^®’^^ that follows the final adiabatic stage 
and that has no driving at all. 

In the following, we shall first present the scaling the¬ 
ory and obtain different scaling behaviors in different 
stages and their crossovers. These are then confirmed 
by simulations on a two-dimensional (2D) classical Ising 
model. As this is a generic model for critical phenomena, 
we expect the results to be applicable to other models and 
even to quantum critical behavior as well. We shall not 
study the topological defects as their counting is not easy 
and detecting scaling behavior of other observables has 
been advocated . 

Near the critical point, the scaling behaviors of macro¬ 
scopic quantities can be readily described by a scale 
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FIG. 1. (Color online) Evolution of the order parameter M 
for the two-dimensional classical Ising model. The initial dis¬ 
tance to the critical point eo and the initial order parameter 
Mo are indicated except the one with eo = —1, where any 
given Mo will always decay rapidly to the equilibrium values. 
Heating instead of cooling is performed to facilitate presenta¬ 
tion. Different stages are marked. 

transformation. Our scaling theory is based on 

M{t,R,Mo,eo,e) 

= C/(Mo, 6), 

for a rescaling of a factor b, where |£o| -C 1, /3 is the crit¬ 
ical exponent for M, and U{Mo,b) is the universal char¬ 
acteristic function describing the rescaled initial magne¬ 
tization®^’®'^. In Eq. (I), we have purposely written t, 
i?. So, and e out though they are not independent as 
£ = £0-1- Rt. For a small Mq, U{Mq, b) = Mob^° with xq 
being the scaling dimension of U{Mo, b) has one 

fixed point U{0,b) = 0 for arbitrary b > 1. For a hard- 
spin system, in which M is bounded, the saturated Mq 
is another fixed point, since the rescaled Mq is invariant 
under coarse graining®®. Note that the two additional 
scaling variables, £o and Mq, are present only for small 
|£o|; for large |£o|, they are absent as they are then irrel¬ 
evant. Equation (1) with a small Mq can be justified by 
an RG theory which combines the critical initial slip®®’®^ 

and the FTS theory®^, 57,58^ 

The scaling forms of different stages can now be ob¬ 
tained from Eq. (1) by comparing the relevant time 
scales. In the first stage, in which t is small, C,i is small 
and growing. Accordingly, relaxation dominates. By set¬ 
ting tb~^ = 1, we arrive at the scaling form 

M(t, R, Mo, £o) = C/(Mo, 

( 2 ) 

where fi is a scaling function. It is valid when all scaled 
variables are small. In particular, RV^/^ <C 1, or t <C 
R-^/^'^, i.e., (^i Cd as ought to be. Detailed scaling 
behavior can be obtained from Eq. (2) as follows. 

Eor £o = 0 and a small Mq, 

One can expand fi in RV^^^ and to the second 

order and obtains 

M ~ Mot«A(0,0,0) (0,0,0). (3) 

where 9 = (xq — (3jv)!z and a prime stands for a par¬ 
tial derivative. In Eq. (3), the first term describes the 
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usual critical initial slip''’^’®'*^. The second term of Eq. (3) 
displays the driving-induced deviation from the critical 
initial slip. It is a mixed term between Mg and R and 
arises from the fact that, if Mg = 0, M remains zero as 
e does not break the symmetry. The external driving, 
which dominates near the critical point in the ordinary 
KZM, here acts only as a perturbation. 

For Eg = 0 and the saturated Mg, U{Mo,t^^^) = Mg. 
In the initial stage, M now decays according to 

M ~ t-^/'^Vi(0,Mg,0) -bi?t('^’'^-^)/'^Vi(0,Mg,0), (4) 

where the first term is the nonequilibrium relaxation®^ 
and the second term arises again from the perturbation 
of the driving. Note that for i? > 0, /i(0,Mg,0) < 0 
because M must decrease as the temperature increases. 

Crossover to the FTS stage occurs at ^ 1, or 

ii ~ ^ (d- This is not i of the KZM as it is the 

crossover from the relaxation stage, which is also nonadi- 
abatic. However, the asymptotically identical forms show 
that (d of the FTS regime does not depend on the initial 
conditions. The scaling form of the FTS stage can be 
obtained from Eq. (1) as 

( 5 ) 

with another scaling function / 2 . For eg = 0 and the 
saturated Mg, Eq. (5) is quite similar to the usual FTS 
form®®’®^’®®. However, the scaling functions are differ¬ 
ent, because they characterize different evolutions from 
distinct initial conditions as can be seen from Fig. 1. 

When » 1, or "C Cd, the system enters 

the adiabatic stage with a scaling form 

M(i?,Mg,eg,e) = £^/ 3 (eg£“\C/(Mg,e“''),i?e“‘'’"^), 

( 6 ) 

where /g is a scaling function. This crossover is similar 
to the usual impulse-adiabatic crossover in KZM as can 
be seen in Fig. 1. Indeed, the time when the curve of the 
usual KZM tends to zero is close to the corresponding 
time of the curve starting with an nonequilibrium state. 
This indicates again that the timescale in the FTS stage 
is consistent with the timescale in the impulse region. 

The scaling theory is applicable to other situations in 
which other variables than T are changed starting with 
a nonequilibrium initial state near the critical point. For 
example, consider changing the symmetry-breaking field 
h as h = ho + R^t with a small hg and a constant Rh- 
We set £ = 0 to reduced competing scales. In this case, 
there exist also three stages in the driving process. In 
the relaxation stage, the scaling form for small Mg is 

M = Mg^"«/^ (7) 

with Th = z + while in the FTS stage in which 

Ci ^ Cd ~ the scaling form changes to 

( 8 ) 




FIG. 2. (Color online) Three stages of the evolution 
of M under increasing e with fixed = 0.2 and 

= 0.06 for three R indicated. The curves be¬ 
fore and after rescaled are shown in (a) and (b) respectively. 
Semi-logarithmic scales are used. 


where fih and f 2 h are scaling functions. Finally comes 
the h dominated adiabatic stage. Again, f 2 h for /ig = 0 
and a saturated Mg is different from the usual one with 
adiabatic initial conditions. We note that this case has 
been considered in Ref. 42, where a method to deter¬ 
mine the critical exponents was proposed. However, the 
relaxation has not been discussed there. 

To confirm the scaling theory, we take the 2D classical 
Ising model as an example. Its Hamiltonian is 

H = - (9) 

<ij> i 

where 5^ = ±1 and the first sum is over all nearest 
neighbors and the second over all spins. Note that un¬ 
less changing the symmetry-breaking external field h, 
we set h = 0 for simplicity. The critical point of (9) 
is Tc = 2/log(-\/2 + 1)"'^ and the critical exponents 
are (3 = 1/8, z/ = 1, 5 = 15“^^ 2 = 2.1667®®, and 
9 = 0.191®®“®®. They will be taken as inputs to verify the 
scaling forms. The single-spin Metropolis algorithm®® is 
used. The lattice size is 5000, which has been checked to 
produce negligible size effects. Periodic boundary con¬ 
ditions are applied throughout. We calculated averages 
over between 2000 and 3000 samples, which guarantee 
that the relative uncertainty is smaller than 1%. The 
initial conhguration is a uniformly-distributed random 
assignment of S'i = ±1 with an average equal to Mg. 

Firstly, we classify the different stages of the evolution 
and examining the scaling form (5) for small Mg. Fig¬ 
ure 2 shows the dependence of M on e for several Mg 
and £g > 0. When e is small, M increases with e and 
thus t at short times. This is similar to the critical ini¬ 
tial slip in the pure relaxation and is thus the relaxation 
stage. When e gets larger, M decreases as e increases. 
Yet, M increases with R and hysteresis occurs. This is 
the generic behavior of FTS stage®^’®®. Then follows the 
adiabatic stage, in which M is zero, independent of R 
and the initial condition. Because Mg and £g have been 
chosen in such a way that Mgi?“^°/’’'^ and are 

fixed, the curves collapse onto each other after rescaling 
according to Eq. (5), confirming that Mg and eg are in¬ 
dispensable scaling variables. 
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FIG. 3. (Color online) Difference of M between R ^ 0 and 
R = 0 with (a) a small Mo and (b) the saturated Mo for 
So = 0. Double-logarithm scales are used. The insets plot 
the original curves, from which one sees that the initial slip 
emerges after ten Monte Carlo steps per spin or so. 


Secondly, we study the effects of the external driving 
according to Eqs. (3) and (4). In Fig. 3(a), the difference 
between the driving relaxation and pure relaxation satis¬ 
fies a power-law relation, |M —(0, 0,0)| oc 
according to Eq. (3). The fitted slope is 6 + rrjz = 
1.696(2), which agrees with the theoretical value of 
d + "Tt/z = 1.652. For the case of the saturated Mg, 
which is Mg = 1 for Ising model, as shown in Fig. 3(b), 
\M — fi(0, 1, 0)1 changes with t with an exponent 

1.499(4), which is close to {wt — P)/>zz = 1.403, consis¬ 
tent with Eq. (4). The deviations arise from the contri¬ 
butions of higher order terms in the expansions. 

Thirdly, we further verify the scaling theory by exam¬ 
ining the scale transformation (1) for large Mg. In this 
case, the rescaled initial order parameter Mq = [/(Mg, b) 
is not a simple power-law®^’®"^. So, for a given b and 
Eg = 0 for instance, we first estimate Mq from the pure 
relaxation by select an Mg starting with which the evo¬ 
lution of M matches, after its M and t being rescaled by 
b~^/’^ and 6“^, respectively, that starting with Mg®^’®'^. 
With this Mg, the evolution of M when e is changing 
again matches well that starting with Mg upon proper 
rescaling, including R = Rb'"^, as is illustrated in Fig. 4. 
Also manifest in the figure is the three stages similar to 



FIG. 4. (Color online) The evolution of M under changing 
e with £0 = 0 and Mq = 0.2. It matches the rescaled curves 
with the estimated value of Mq and R — Rlf^ for b = 4. 
Semi-logarithmic scales are used. 
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FIG. 5. (Color online) (a) Three stages of the evolution 
of M under changing h with fixed = 0.01 and 

MqR'^^°^^’' = 0.04 for three Rh indicated, (b) The rescaled 
curves. Semi-logarithmic scales are used. 


the case of small Mg. These show that the effects of driv¬ 
ing and of the initial conditions are independent and thus 
confirm Eq. (1). 

Fourthly, we consider the situation of changing the 
symmetry-breaking field h. Figure 5 shows the results 
of changing h as h = ho — R^t with some small and neg¬ 
ative ho. The three stages also show manifestly similar 
to those in Fig. 2. The rescaled curves with different Rh 
and Mg collapse well onto each other for fixed 
and This confirms both that the scaling 

form must include ho and Rh as scaling variables as 
Eq. (8) indicates and that the relaxation-FTS-adiabatic 
scenario is generally applicable in the driving dynamics 
with nonequilibrium initial states near the critical point. 

In summary, we have systematically studied the driv¬ 
ing dynamics starting with a nonequilibrium initial state 
near the critical point. This initial condition breaks 
the adiabaticity and thus changes the adiabatic-impulse- 
adiabatic scenario of the KZM into the relaxation-FTS- 
adiabatic scenario by suppressing the initial adiabatic 
stage. A scaling theory that combines FTS with critical 
initial slip has been developed and account well for the 
universal scaling behavior in this nonequilibrium nonadi- 
abatic case. Numerical simulations on the 2D Ising model 
have confirmed that the theory applies well both to vary¬ 
ing temperature and to varying the symmetry-breaking 
external field. Our theory might provide a way of nona- 
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diabatic quantum computations as opposed to the adi¬ 
abatic ones'*^, as one may now quench nonadiabatically 


from the ground state of an initial Hamiltonian to the 
targeted one even at the critical point of the latter. 
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